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Abstract 


Hardfill dams are comprised of cemented granular materials gathered 
directly from the river bed or its surroundings, or which are by-product 
of the dam excavation. This fact leads to significant variations in as-built 


strengths and deformability parameters within the dam body. In this 
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paper, a numerical approach is applied to derive fragility curves for 
hardfill dams under seismic loading. Fragility curves describe the 
probability of a structure to experience a certain damage level for a 
given earthquake intensity measure, providing a relationship between 
seismic hazard and vulnerability. Both static and seismic responses of 
the hardfill dam are assessed using tridimensional finite difference 
analyses carried out with the software FLAC3D, considering increasing 
levels of seismic intensity. The fragility curves are estimated in terms of 
peak ground acceleration at the free field, based on the evolution of 
damage with increasing earthquake intensity. The proposed fragility 
models allow to characterize the seismic risk of representative hardfill 
dam typologies and foundation conditions considering the associated 
uncertainties, and partially fill the gap of data required in performing a 


risk analysis of this type of dams. 


Keywords: Hardfill dams, fragility curves, vulnerability, risk analysis. 


Resumen 


Las presas de hardfill están construidas con materiales granulares 
cementados recolectados directamente del cauce del río y sus 
alrededores, o producto de excavaciones para la construcción de la 
presa. Este hecho conduce a variaciones significativas en la resistencia y 
los parámetros de deformación en la construcción del cuerpo de la 


presa. En este artículo se aplica un enfoque numérico para obtener 
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curvas de fragilidad de presas de hardfill bajo carga sísmica. Las curvas 
de fragilidad describen la probabilidad de una estructura a experimentar 
un cierto nivel de daño para una intensidad de sismo dada, 
proporcionando una relación entre el riesgo sísmico y la vulnerabilidad. 
Las respuestas estática y sísmica de la presa de hardfill se evalúan 
usando un análisis tridimensional de diferencias finitas realizado con el 
software FLAC??, considerando diferentes niveles crecientes de 
intensidad sísmica. Las curvas de fragilidad se estiman en términos de la 
aceleración máxima del terreno en el campo libre, PGAfsr, basadas en la 
evolución del daño con el incremento en la intensidad sísmica. Los 
modelos de fragilidad propuestos permiten caracterizar el riesgo sísmico 
de presas de hardfill con esta tipología y de las condiciones del suelo, 
considerando las incertidumbres asociadas, y llena parcialmente la 
brecha de datos requerida para realizar un análisis de riesgo de este tipo 


de presas 


Palabras clave: presas de hardfill, curvas de fragilidad, vulnerabilidad, 


análisis de riesgo. 
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Introduction 


Concrete-faced hardfill dams are built mostly using granular materials, 
such as sand and gravel, obtained directly from the riverbed or its 
surroundings, or the same dam excavation by-products, artificially 
cemented together. It is common to have hardfills with 5 to 10% 
content of non-plastic fines. In recent years, the use of hardfill dams 
have increased worldwide due in part to its low cost and simplicity of 
construction. Generally, a typical hardfill dam is constructed using 
symmetrical trapezoidal sections. An impervious concrete-face is placed 
on the upstream slope. Hardfill dams heights range from 3 to 110 m 
approximately (Cai, Wu, Guo, € Ming, 2012; Mason, Hughes, 8 
Molyneux 2012; Fujisawa € Sasaki, 2012; Guillemot € Lino, 2012). 
Traditionally, they have been classified as small to medium size dams 
(Cai et al., 2012; Guillemot € Lino, 2012). In particular, one of the 
highest hardfill dams reported in the literature is the Cindere dam, 
located in Turkey, which is 107 m high, and has a capacity of 84.3 
million m* (Guillemot 8£ Lino, 2012; Batmaz, 2003). The trapezoidal 
shape helps minimizing stresses within the dam body and foundation as 
well as reducing oscillations during a seismic event (Xiong, Weng, € He, 
2013). Some of the main advantages of hardfill dams are the quick and 


simple construction because of the strength demand on foundation is 
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relatively less than that required in other types of dams such as gravity 
dams with typical sections, and also that the specifications of 
aggregates and cement are less stringent, which results in reduced 
production cost of hardfill. These facts, however, can lead to significant 
variations in as-built strengths and deformability parameters within the 
dam body (Batmaz € Gúrdil, 2003). 


Fragility curves relate the seismic intensity parameter to the 
probability of reaching or exceeding a given level of damage for each 
element at risk. The level of shaking can be quantified using numerous 
earthquake parameters, including peak ground acceleration (PGA), peak 
ground velocity (PGV), permanent ground deformation (PGD), spectral 
acceleration, spectral velocity or spectral displacement. To build these 
curves it is necessary to define some objective indicators that allow to 
estimate when the damage of a structure has passed from one level of 
damage to another. These levels are known as damage thresholds or 
damage degrees (e.g. minor, moderate, extensive and complete). 
Fragility functions are required elements for the assessment of 


vulnerability and risk (Argyroudis 8 Kaynia, 2015). 


Fragility curves are essentially based on expert opinions and field 
observations. Nevertheless, one  limitation when  dealing with 
geotechnical problems, is the lack of well documented field data, with 
respect to damage indicators (e.g. permanent displacement, cracks and 
resisting exceedance force). The objective of this paper is to develop 


numerically-derived fragility curves for hardfill dams following the 
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procedure proposed by (Argyroudis € Kaynia, 2015). Thus, commonly 
used geometry, structural parameters, rock foundation conditions, and a 
wide variety of ground motion characteristics are considered. This 
information partially fills the gap of information regarding the seismic 


response of this type of dams under medium to strong ground shaking. 


Methodology 


The procedure for developing numerically-derived fragility curves is 
described in Figure 1. First, the dam typology and hardfill and rock 
foundation mechanical parameters must be established, along with the 
material properties of the concrete face. Second, the damage states, ds, 
are stated, and the dynamic analysis of the hardfill-rock foundation dam 
is performed, considering a three-dimensional seismic environment. 
Third, the damage evolution is obtained for each intensity measure, MI, 
defining the value of the average threshold of MI for each, ds. Finally, 
fragility curves are derived for four damage states (i.e. minor, 
moderate, extensive and complete). These curves consider the 


uncertainty in seismic demand, capacity of each element, and the 
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definition of the damage index and damage states. The response of the 
basin-rock foundation-dam coupled system is computed through a series 
of 3D finite difference dynamic analyses for an increasing level of 
seismic intensity using the software FLAC*” (Itasca Consulting Group, 
2005). The geometry, material properties, and structure details are 
parameters, which describe the typology of the haradfill dam and its 
capacity to withstand seismic loads. The seismic loads are function of 
the seismic environment and basin-rock foundation-dam interaction. 
Three types of rock foundation were considered in the numerical study, 
to assess the effect of different deformability parameters on the 
vulnerability of a harafill dam. Several ground motion were selected in 
terms of amplitude, frequency content and duration. Damage was 
quantified in terms of the shear stress acting on the dam under static 
and seismic loading, normalized by the shear strength (stresses 


capacity) on the dam. 
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Hardfill dams typology Rock foundation 


Geometry, materials and design Typical strength and 
parameters deformability parameters 


Damage index (DI), Seismic input motion 
damage states (ds), 3D dynamic analysis of Accelerograms, intensity levels 
thresholds values of DI hardfill dams 
for each -ds 


Uncertainties 
(Seismic demand, element 
capacity, definition of Dl and ds) 


Evolution of damage with intensity measure 
(IM), definition of median threshold value 
of IM for each ds 


Fragility curves for a hardfill dam 
with different types of rock foundation 


Figure 1. Procedure followed herein for deriving numerical fragility curves for harafill 


dams. 


Numerical model 
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A numerical study was conducted in order to characterize the damage 
that a typical hardfill dam can exhibit for several seismic shaking 
scenarios varying from moderate to extreme ground shaking. It is worth 
mentioning that all the records used correspond to earthquakes 
measured at rock outcrops. A schematic representation of the dam 
geometry is presented in Figure 2. The dam is symmetrical, and 107 m 
high. The faces slopes are 0.7:1.0. The dam is considered founded on 
rock. Table 1 summarizes information regarding the earthquakes 
considered in the analyses, and Figure 3 shows the corresponding 


response spectra. 
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Figure 2. Dam typology a) Cross section b) Downstream view (dimensions in meters). 
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Figure 3. Response spectra of ground motions considered to derived the fragility 


curves. 
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Table 1. Earthquakes considered in the analyses. 
z z Moment  Epicentral 
Seismogenic Earthquake j Ñ PGA | Duration 
ñ Year |magnitude, distance, (9) ( ) 
ame sec 
Zone M,, Ed (km) g 
Montenegro 
(former 1979 6.9 21 1.774 48 
Yugoslavia) 
Kocaeli (Gebze, 
Normal 1999 7.4 47 0.218 32 
Turkey) 
Umbria Marche 
(Gubbio-Piana, 1998 4.8 10 0.235 40 
Italy) 
Maule (Concepcion 
2010 8.8 109 0.605 120 
San Pedro, Chile) 
Honshu (Haga, 
Subduction 2011 9.0 283 0.814 300 
Japan) 
Michoacan 1985 8.1 265 0.032 170 


(Tacubaya, Mexico) 


Three sets of rock foundation properties were considered in the 


analyses, as compiled in Table 2. The concrete deformation parameters 


used in the analyses are summarized in Table 3. For a given project, it is 


necessary to perform both laboratory and field tests to characterize 


strength and deformability parameters of the hardfill for short and long- 
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term conditions, to conduct a final dam reliability assessment. For the 
cases studied included herein, these properties were defined from a 
statistical analysis of the information available in the technical literature 
(Zou, Li, Xu, 2011; 2008) as 
summarized in Table 4. 


8 Kong, Xiong, Young, € Peng, 


Table 2. Properties of the rock foundation used for the analysis. 


Properties Units Type C1 Type C2 Type C3 
Young modulus MPa 4903 9806 15475 
Shear Modulus, G MPa 2114 4227.0 6670 
Volumetric weight y kN/m* 25 26 27 
Poisson ratio v 0.16 0.16 0.16 


Table 3. Concrete properties used in the analyses. 


Properties Units Concrete 
Young modulus, E MPa 28000 
Volumetric weight, y kN/m* 24 
Poisson ratio, v 0.20 


Table 4. Typical hardfill properties reported in the technical literature. 


] Friction Young . o. 
Cohesion, Cement Curing | Confining 
angle, 4 | modulus, E . Reference 
c (MPa) content time stress 
(9) (MPa) 
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(kg/m*”) | (days) (kPa) 
0.80 49 342 198 14 100 
0.42 46 341 132 14 100 Zou et al., 2011 
0.25 44 145 66 14 100 
0.24 49 594 60 28 200 
0.20 49 500 50 28 200 Sun 8 Yang, 
0.12 49 450 40 28 200 2011 
0.08 46 400 30 28 200 
0.08 42 250 20 28 200 
0.5 45 2,000 60 91 =-- Xiong et al., 2008 
4,000- Yanmaz €: Sezgin, 
0.8 45 70 180 -- 
11,586 2009 


Short-term condition analyses 


For short-term conditions, the hardfill cohesion, c, corresponds to the 
average of those values associated to a curing time of 14 days gathered 
from Zou, Li, Xu and Kong (Zou et al., 2011; Xiong et al., 2008) 
summarized in Table 4. Similarly, the internal friction angle, oe, 


corresponds to the lowest reported value. 
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Long-term conditions analyses 


For the long-term conditions, the c and q values associated to a hardfil! 
curing time of 180 days were considered (Sun € Yang, 2011), which 
correspond to the mechanical properties of the hardfill parameters used 
in Cindere dam (Table 4). To take into account the evolution of hardfil! 
properties with curing time a parabolic distribution was assumed. Figure 
4 shows the evolution of the hardfill friction angle, cohesion and initial 
tangent elastic modulus considered in this study. The short and long 


term conditions parameters are summarized in Table 5. 
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Figure 4. Evolution a) Friction angle, b) Cohesion and c) Initial tangent elastic 
modulus of harafill with curing time. 


Table 5. Short and long term conditions parameters. 


Condition | Cohesion, | Friction angle, Initial Volumetric | Poisson 
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c (MPa) 9 (9) Young weight, y ratio, v 
modulus, (kN/m?) 
Eo (MPa) 
Short-term 0.5 42 420 23 0.20 
Long-term 0.8 45 7930 23 0.20 


Tridimensional soil-structure interaction 


analyses 


The seismic dam response analyses were carried out with a three- 
dimensional finite difference model, using the software FLAC*” (Itasca 
Consulting Group, 2005). The numerical model is depicted in Figure 5. 
The dam-foundation-embankment system was modelled with solid 
elements. This allows obtaining acceleration, velocity and displacements 
time histories simultaneously, along with strains and stresses. An 
elasto-plastic Mohr—-Coulomb model was used to represent the stress- 
strain relationship for the hardfill. The rock and hardfill dynamic shear 
stiffness and damping variation as a function of shear strain during the 


seismic event was assumed negligible, considering the soundness of the 
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rock foundation and that the hardfill should behave in its elastic range to 
avoid excessive and unacceptable cracking. The calculation is based on 
the explicit finite difference scheme, to solve the full equations of 
motion, using lumped grid point masses derived from the real density of 
surrounding zones. The added masses, required to account for 
hydrodynamic pressures in the numerical simulation, were obtained 
using the equation proposed by Chwang and Housner (1990). In the 
three-dimensional model a typical hardfill construction procedure was 
considered, which consists in a continued placement of hardfill layers 
with constant thickness (i.e. 1 m). The placement of one haradfill layer 
per day was simulated, considering the evolution of the hardfill 
properties with curing time explicitly in the model. The dam stages 
modelled were: 1) Sequential dam construction, 2) First reservoir filling, 
and 3) Seismic loading. The model depth is 214 m, the width is 480 m 
and the length 1080 m. The free field boundaries implemented in FLAC?? 
(Itasca Consulting Group, 2005) were used through the edges of the 
model. A flexible base was considered along the bottom of the model, to 
avoid boundary seismic wave reflection towards the  dam-rock 


foundation system. 
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Boundaries 


Rock 
Boundaries 


214 m 


1080 m 


478.4 m 


Figure 5. Finite difference model (FLAC*D). 


Fragility curves development 


Hardfill dam numerically-derived fragility curves were obtained following 
the approach described in Figure 1. As previously described, the dam 


was assumed founded on rock. Three cases of rock foundation were 
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considered to take into account the influence of their deformability 
parameters. The seismic loads are function of the surrounding rock, 
seismic environment and the basin-rock foundation—-dam interaction. 
Several ground motions were selected based on amplitude, frequency 
content and duration. Finally, fragility curves were derived for different 
damage states (minor, moderate, extensive and complete) considering 


the primary sources of uncertainties. 


Both normal and subduction earthquake type events were 
considered. A total of six ground motions recorded during these events 
were used in the analyses, exhibiting varying spectral acceleration 
amplitudes, frequency contents, significant durations and seismotectonic 
settings (Figure 3). The time histories are scaled to five intensity levels 
of PGA (i.e. 0.15, 0.30, 0.45, 0.60 and 0.759) in order to perform the 
dynamic analyses of the basin-rock foundation-soil system for an 


increasing seismic excitation. Thus, a total of 90 cases were analyzed. 


Damage index and damage states 


A practice-oriented damage criterion for hardfill dams was proposed. 


Considering that there is currently a lack of relevant information on this 
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subject, in this research the damage index, DI, is defined as the 
relationship between the shear stress acting on the dam, Tact, 
considering both static and seismic loading conditions, normalized by 


the shear strength (stresses capacity), Tcap, as follows: 


DI = “22 (1) 


Tcap 


where: 


== y(01-02)2+(0,-03)2+(03-01)2 
o al E 


(2) 


Tcap = C+ (Oj. tano) (3) 


c is the hardfill cohesion, q is the friction angle, the octahedral stress, 


(0,+0,+03) 


Ooci 1s defined as , and 0,,0,,03 are the principal stresses. 


Four damage states were considered (i.e. minimum, moderate, 
extensive and complete, as summarizes in Table 6. In particular, Figure 
6 shows the damage index contour through the dam body for the 
seismic subduction event of Maule, Chile, 2010, scaled to 0.75 g and 


considering the rock type 1. It can be observed that, in this case, the 
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higher values of damage indices are concentrated in the downstream 


part of the dam. 


Table 6. Definition of damage states for the hardfill dam. 


Damage states 
(dsi) 


dsi 


Minimum 


ds2 


Moderate 


ds3 


Extensive 


ds4 


Complete 


Range of damage 


index 


0.5< DI 
<0.65 


0.65< DI 
<0.8 


0.8< DI 
<1.3 


1.3>DI 


(a) 
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Figure 6. Damage Index contour (a) downstream face and (b) upstream face for the 


seismic subduction event of Maule, Chile, 2010. 


Derivation of fragility curves 


Fragility functions describe the probability of reaching or exceeding a 
given damage state for a certain intensity measure. They are usually 


described by a lognormal probability distribution function, as follows: 
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_ 1 PGAsf 
PA(PGAff) = 0 (on (Az) (4) 


where Pf e) is the probability of exceeding a particular damage state, ds, 
for a given seismic intensity level defined by the earthquake intensity 
measure, IM (e.g.; peak ground acceleration, PGA), O is the standard 
cumulative normal distribution, PGAfsfmi is the median threshold value of 
the earthquake intensity measure, required to cause the ith damage 
state and Btot is the total lognormal standard deviation. Therefore, the 
development of fragility curves requires the definition of two 
parameters, PGAfsfmi and Bro (Mayoral, Argyroudis, € Castañon, 2016). 
The lognormal standard deviation Bro, which describes the total 
variability associated with each fragility curve, is estimated considering 
the uncertainty in the definition of damage states (Bas), the response 
and resistance (capacity) of the element (fc), and in the earthquake 
input motion (demand) (Ba). Due to the lack of a rigorous estimation, 
Bas is set equal to 0.4 following HAZUS (NIBS, 2004) approach for 
buildings, Bc is considered equal to 0.4 based on engineering judgment 
(Argyroudis 8 Kaynia, 2015), while fB¿ is estimated based on the 
variability in the hardfill response (DI) that have been calculated for the 
selected ground motions. The value of Btot is estimated as the root of the 
sum of the squares of the component dispersions. The median threshold 


value of the earthquake parameter, PGAfrmi is obtained for each damage 
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state based on the evolution of damage with increasing earthquake 


intensity. 


Regression analysis 


A regression analysis is performed, considering the damage index (DI) 
as the dependent variable, and the peak ground acceleration at free 
field, PGAff, as the independent variable, as it is illustrated in Figure 7, 
Figure 8 and Figure 9. The estimated parameters of the fragility curves 
are given in Table 7, and the corresponding fragility curves are 
illustrated in Figure 10, Figure 11 and Figure 12. The fragility curves 
were developed based on the peak ground acceleration determined at 
the free field, PGAfsf, for both subduction and normal-fault events scaled 
up to 0.75 g of PGAff. Evidently, the derived fragility curves can be 
applied for preliminary purposes only to cases where the dam 
foundation, foundation rock conditions and seismic environment are 
similar to those considered in the analyses. The development of this 
kind of functions, is a key step in seismic risk analyses to capture the 


most important seismic failure modes of the dam. Also, these are 


Tecnología y ciencias del agua, ISSN 2007-2422, 11(1), 132-168. DOI: 10.24850/j-tyca-2020-01-04 


157 


Tecnología y 


CienciaszAgua 


2020, Instituto Mexicano de Tecnología del Agua 


Open Access bajo la licencia CC BY-NC-SA 4.0 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


required to perform a proper seismic vulnerability assessment to 


foreseen the capacity of the system to withstand very large to extreme 


events, ensuring earthquake preparedness, and reducing life losses or 


post-earthquake distress. 


Table 7. Parameters of the fragility curves. 


Case 1 


Case 2 


Case 3 


Median threshold 


Median threshold 


Median threshold 


PGA ffmi Bro PGA ffmi Bro PGA ftmi Brot 
(8) (8) (8) 

Minor 0.47 0.569 0.76 0.614 1.10 0.572 
Moderate 2.51 0.569 3.36 0.614 6.62 0.572 
Extensive 5.01 0.569 6.54 0.614 13.39 0.572 
Complete 8.24 0.569 10.63 0.614 22.10 0.572 
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DI type soil C1 
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Figure 7. Evolution of DI with PGA for rock type C1. 
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Figure 8. Evolution of DI with PGA for rock type C2. 
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Figure 9. Evolution of DI with PGA for rock type C3. 
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Figure 10. Fragility curves for Hardfill dam for rock type Cl. 
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Figure 11. Fragility curves for Hardfill dam for rock type C2. 
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Figure 12. Fragility curves for Hardfill dam for rock type C3. 
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Probability of failure 


Uniform hazard spectra for four return periods (i.e. 125, 250, 475 and 
2475 years) were obtained from the seismic hazard curves for both 
normal and subduction fault events. As it is well known, the uniform 
hazard spectra, UHS, is a representation of the relationship between the 
natural vibration period, T, and spectral acceleration, Sa, for a given 
exceedance probability associated with a return period. The uniform 
hazard spectra are shown in Figure 13a and Figure 13b for normal and 
subduction events respectively, which were used to obtain the peak 
ground accelerations (PGA). Then, with the PGA, and the numerically- 
derived fragility curves, the probability of reaching or exceeding the 
different states of damage (i.e. minor, moderate, extensive and 
complete) was computed, for each rock foundation type C1, C2 and C3, 
and the results are shown in Table 8. Overall, the numerically-derived 
fragility curves lead to high values of probability of damage. This 
appears to be over conservative, thus, more research is needed on the 


definition of the damage states to ensure more realistic results. 
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Figure 13. Uniform hazard spectra for Tr = 125, 250, 475 and 2475 years (a) normal 


Table 8. Probability of damage for C1, C2 and C3. 


and (b) subduction events. 


Probability of Damage (%) 


Seismogenic | Period | PGAsock 
ci Cc2 C3 
Zone Return (9) 
Mi? Mo? E? c* Mi: | Mo? E? c* Mi: | Mo? E? C* 
125 5 
years : 0.391 | 0.159 | 0.056 | 0.008 | 0.939 | 0.525 | 0.194 | 0.025 | 0.714 | 0.164 | 0.028 | 0.001 
250 $e 
Normal years " 0.583 | 0.304 | 0.136 | 0.028 | 0.980 | 0.719 | 0.365 | 0.074 | 0.857 | 0.317 | 0.080 | 0.006 
475 07 
years " 0.623 | 0.341 | 0.159 | 0.035 | 0.985 | 0.755 | 0.407 | 0.090 | 0.879 | 0.355 | 0.097 | 0.008 
2475 0.22 0.773 | 0.511 | 0.287 | 0.084 | 0.996 | 0.876 | 0.590 | 0.191 | 0.948 | 0.531 | 0.197 | 0.026 
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years 
125 mm 
years : 0.335 | 0.126 | 0.041 | 0.005 | 0.918 | 0.462 | 0.154 | 0.017 | 0.660 | 0.129 | 0.020 | 0.001 
250 yA 
years : 0.391 | 0.159 | 0.056 | 0.008 | 0.939 | 0.525 | 0.194 | 0.025 | 0.714 | 0.164 | 0.028 | 0.001 
Subduction 
475 e 
years j 0.493 | 0.230 | 0.092 | 0.016 | 0.966 | 0.633 | 0.279 | 0.045 | 0.798 | 0.239 | 0.051 | 0.003 
2475 7 
years : 0.623 | 0.341 | 0.159 | 0.035 | 0.985 | 0.755 | 0.407 | 0.090 | 0.879 | 0.355 | 0.097 | 0.008 


Conclusions 


A numerical approach is applied to construct fragility curves for hardfill 
dams under seismic loading. The response of the hardfill dam is 
evaluated based on 3D numerical models using the software FLAC??, 
Seismic inputs with several frequency content and scaled to various 
levels of seismic intensity (PGAfFÉf) and typical rock profiles (Table 2) are 
considered. A practice-oriented damage criterion for hardfill dams was 
proposed. The damage index (DI) was established as the ratio between 
the shear stress acting, considering both static and seismic loading 
conditions, normalized by the shear strength. The fragility curves were 


derived as a function of the peak ground acceleration in the free field, 
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considering the related uncertainties in the definition of damage states, 
the demand, and the capacity of the hardfill material used to build the 
dam. The applicability of the derived fragility curves has some 
constraints, which are related to the specific assumptions of the 
typology of the dam and selected parameters in the numerical models, 
so they become very specific for the cases studied, having to develop 
fragility curves for each condition and dam typology. Therefore, the 
provided fragility curves can be applied for preliminary purposes only, in 
those cases where the structure and rock conditions are similar to those 
considered in the analyses. Overall, the numerically-derived fragility 
curves lead to high values of probability of damage. This appears to be 
over conservative, and more research is needed on the definition of the 


damage states to ensure more realistic results. 
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